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ABS n T developments in two-equation models and second order closure models wUl 
be presented Two-equation models (e.g., k-e model) have been wrdely used m CFD 
engineering problems. Most of low-Reynolds number two-equation modds contarn some 
wall-distance damping functions to acount for the effect of wall on tur u ence^ 
this often causes the confurions and difficulties in computing flows wrth comp ex geomet^ 
and also needs an ad hoc treatment near the separation and reattachment pom • 
this paper, a set of modified two-equation models is proposed to remove abovement.oned 
swcomi^. The calculations using various two-equation models are compared wrth 

direct numerical simulations of channel flows and flat boundary layers. 

Development of second order closure model will be also discussed with emphases on the 
modeling of pressure related correlation terms and dissipation rates in the second momen 
equations. All the existing models poorly predict the normal stresses near t e w an 
to predict the 3 dimensional effect of mean flow on the turbulence (e.g decrease 
shear stress caused by the cross flow in the boundary layer). The new y 
order near-wall turbulence model to be described in this paper rs capable of captunng h 
near-wall behavior of turbulence as well as the effect of three dimensron mean flow 


turbulence. 
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1. fc-e model 

The two-equation model, especially k-t model, is still the most widely used model 
for computing engineering flows. We first list some nf i 

models, J d lheir oredii T ” , , * US " 1 

, , ’ „ thClr predlctlons on the f ««y developed channel flows and 

un ary ayer flows compared with the corresponding direct numerical simulations Then 

nZT T1 k ' e m ° del WhiCh dOCS n0t — - distance. The proposed 

model has been also tested using direct numerical simulation data. 

The eddy viscosity v T is assumed in two-equation models as follows: 

k 2 

V T = C M f M — 


or 


V T = Cftf^kr 

where r = k 2 jl 

The general k- e (or k-r) model equations are of the following forms: 


k,t + Ujkj = 

e > t + U i e ,j = 
r .‘ + U i T J = ^ 


(o-* + ") *•*] . + 11 + V T ( Ui,j + U jt i) U itj -e + D 
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This table lists the model terms, damping functions and model constants appeared in 
various two-equation models 
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The following figures show the predictions on fully develop channel flow using various 
two-equation models compared with the direct numerical simulation dataJ 7 ! Ploted quan- 
tities include the mean velocity U, turbulent shear stress (uv), turbulent kinetic energy k 
and dissipation rate EPS t. The open circle represents direct numerical simulation, and 
the solid line represents the model prediction. 
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The figures below show the model predictions on flat plate boundary layer flow. A 
direct numerical simulation of boundary layer flow^ is used for comparison. The skin 
friction coefficient C/ is also included in the comparisons. 
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Overall, Shih’s k-e model gives better prediction in both fully developed channel flows 
and flat plate boundary layer flows according to the comparisons with corresponding direct 
numerical simulations. However, this model has the same problem as the others, that is it 
contains the wall distance parameter y+, which is defined as 

+ u ry 

y = — 

V 

where u T is the friction velocity. The difficulty would occure in some situations. For 
example, near the seperation point u T aproaches zero and hence v t (through /„(y+)) will 
approach zero everywhere when this « r is used. Another example is the flow with complex 
geometry that the wall distance is not well defined. In the both cases, the ad hoc treatment 
is needed in the model implimentation. We notice that Jonse-Launder’s modelW does not 
contain y+ However, its present form does not perform very well in the simple testing 
flows. Here, we based on the Shih’s k-e model modify the parts which is related the y+ 
with another parameter i? = R is a ratio of turbulence length scale to viscous 

length scale. The modifications* 9 ) made here are: 

U = 1 ~ exp jc 3 [l - exp(C' 6 f2 1 / 4 )J J 


‘ J iJ 

c = e [l - exp(-i? t 1/2 )] 

where R t = k 2 /ue, C 0 = .004, C 3 = .0004, C 6 = 1.2 
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turbulent viscosity 


The following figures show the predictions from the present model on fully developed 
channel flow compared with other models (including Jonse and Launder’s model W) and di- 
rect numerical simulation dataJ 7 ) The open symbols represent direct numerical simulation, 
and the lines represent the model prediction. 
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2 . Second order modeling of near- wall turbulence 

Using the near-wall asymptotic behavior of turbulence^ 10 ] as model constraints, we 
formed a set of modeled transport equations for the Reynolds-stress tensor and the dissi- 
pation rate of turbulent kinetic energy. The main emphasis was on developing a near-wall 
model for the pressure correlation and dissipation terms in the Reynolds-stress equation. 
A modeled dissipation rate equation is derived more rationally. Asymptotic analysis shows 
that near the wall, the viscous diffusion term in the Reynolds-stress equations becomes 
the leading term and is balanced by the pressure correlation and dissipation terms. We 
use this as a model constraint in the model development. The proposed models satisfy 
realizibility and ensure no unphysical behavior will occur. Here, we briefly describe and 
list the proposed models. 

Reynolds stress equation 

The exact equation for the Reynolds stress tensor is: 

— {tijtij;) = Pij + T{j + D\j ^ + n fi — €{j 

where ( ) stands for an ensemble average, D/Dt = d/dt + U k d/dx *. the terms Pih Tij, 
D\j \ Uij and €ij represent the production, turbulent diffusion, viscous diffusion, velocity 
pressure- gradient correlation and dissipation tensor, and are identified as follows: 


Pij = -{uiU k )U jyk - {ujU^U^k 

Pij — 

= u{uiUj) M 

n,j = + u jp,i) 

£ij = 2i '(u itk u jt k) 

The proposed near-wall model for Ily — €ij is: 


Uij tij 



[2{ti,Uj) + 4({uiU k )rijn k + (ujU^n.n*) + 2(ufcU,)n*n / n;n J ] 


(a 7 ) 2 

where n, is a irnit vector normal to the surface, and f w = exp(“(iJt/Ci) 2 ), Rt = 9 ^ 7*5 
C\ = 1.358i2^ 4 , R^t — u T S/u . u r is the friction velocity, 6 is the thickness of the boundary 
layer or the half width of the channel. 

Away from the wall, the velocity pressure-gradient correlation II , j is split into the 
rapid part Il|^ and the slow part 11 ^: 


n., = n!J» + ng> 
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The proposed model for II - e t j is: 

— £ «i = ~ c (0bij + — fw) 

where 

& = 2 + + 8(U ln t ! + 62.4(— // + 2.3 ///)]} exp(-~) 

F= 1 + 27/// +9// 

n = 4 m J( 

I II = ^ ^ijbjh^ki 

bij = (u lUj )/( 9 2 ) _ ^y/3 

The rapid part of velocity pressure-gradient, n[ 1} is modeled as followsfShih and 
Lumleyt 11,12 !): V 

n l- = ^ + 2a s)(<l 2 )(Ui,} + U jt i) - |(1 - a 5 )(p 0 - |p<!> 0 ) 

+ ^3 + T a ®^' D,i ~ 3 PSij ^ + 15 ^ ~ D 'j} + t bi i p 

2 

+ 5^) {((UiUt) ^ + - ( u i u p){ u j u q)(Up, q + U q ,p)] 

where, 

Pij = ~(uiUk)Uj t t — (ujUk)Ui t k 

Dij = -( UiU k )Uk y j - { UjU k )Uk,i 



a s = -^ (1 + C 2 F 1/2 ) 

C 2 = 0.8[1 - exp(— (Pj/40) 2 )] 

Finally the model for the third moments is modeled as: 

{UiUjUk) = -.°7^[(u t U J ,)(u,U i ) iP + (u iUp )( U ,U t ) p + (UiUpjiujUkj'p] 

Dissipation rate equation 

The modeled dissipation rate equation derived in this work is: 

€, < + Uie ti = (ue ti - (cui)),i - tf’OT^r 

(r) 

- V>i^y - j>2 ^(ukuMUij, - U ltij )Uijk 
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where 


V> 0 = — + 0.98(1 - 0.33 ln(l - 55//)]exp(-2.83i?r 1/2 ) 
5 

xpi « 2.1 
* = -. 15(1 -F) 

. f _Wk 

4{? 2 > 

The turbulent flux term (eu*) is modeled as: 

(cut) = -• 07 ^§7 ( u *«j>) c .J’ 


These figures show some existing Reynolds stress models (for example, Launder and 
ShinuJ 13 !, Lai and So^ 14 ^) and present model compared with the direct numerical simulations.^ 
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3. Second order modeling of a three-dimensional boundary layer 

A study of three-dimensioned effects on turbulent boundary layer were achieved by 
direct numerical simulation of a fully developed turbulent channel flow subjected to trans- 
verse pressure gradient. The results show that, in agreement with experimental datal 16 ^ , the 
Reynolds stresses are reduced with increasing three-dimensionality and that, near the wall, 
a lag develops between the stress and the strain rate. To model these three-dimensional 
effects on the turbulence, we have tried various two equation models and second order 
closure models. None of the current models can predict the reductions in the shear stress 
observed using direct numerical simulations. However, we found that the newly proposed 
second order closure model listed in the previous section do at least qualitatively capture 
these three-dimensional effects. 
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The following figures show the direct numerical simulation of the three dimensional 

boundary layer flow and the model predictions from Launder and Shima, Lai and So and 
the present models. 



Schematic of three-dimensional channel flow. 
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Abstract 

Tins paper presents a set of realizable second order models for boundary free turbulent 
ows. The constraints on second order models based on the realizability principle are re- 
examined. The rapid terms in the pressure correlations for both the Reynolds stress and 
the passive scalar flux equations are constructed to exactly satisfy the joint realizability. 
All other model terms (return- to-isotropy, third moments and terms in the dissipation 
equations) already satisfy realizability (Lumley 1978, Shih and Lumley 1986). To correct 
ie spreading rate of the axisymmetric jet, an extra term is added to the dissipation 
equation which accounts for the effect of moan vortex stretching on dissipation. The test 
flows used m this study are the mixing shear layer, plane jet, axisymmetric jet and plane 
wake The numerical solutions show that the new unified model equations (with unchanged 
model constants) predict all these flows reasonably as the results compare well with the 
measurements. We expect that these model equations would be suitable for more complex 
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Evaluation of Turbulence Models 
for Predicting Buoyant Flows 

Experimental data for the buoyant axisymmetric plume are used to validate certain 
closure hypotheses employed in turbulence model equations for calculating buoyant 
flows. Closure formulations for the turbulent transport of momentum , thermal 
energy, kinetic energy ; and squared temperature used in the k~t and algebraic stress 
models are investigated. Experimental data for the mean velocity , mean temperature , 
and kinetic energy are used in the closure formulation to obtain Reynolds stresses, 
heat fluxes, etc., which are then compared with their measured values. 


1 Introduction 

Various turbulence models have been formulated for pre- 
dicting buoyant-driven flows. Some of the parameters in these 
models have been determined by keying the solution of the 
model equations to experimental data for certain basic flows 
such as decay of grid turbulence. Other parameters have been 
determined by calibrating closure formulations directly with 
experimental data. However, this approach may be somewhat 
inaccurate due to the lack of quality experimental data for 
certain correlations, especially dissipation. Finally, certain 
mode! parameters have been fine tuned or determined by re- 
quiring that the computed solution agree with experimental 
data for more complex flows, such as shear flows. In addition 
there have been instances where model parameters have been 
adjusted or empirical corrective terms added so that agreement 
with experimental data is accomplished for a particular flow. 
When model parameters are adjusted to get agreement, say 
for the mean velocity and temperature fields for a particular 
flow, little regard is given for the internal integrity of the model . 
In other words, are the various processes such as diffusional 
transport, pressure-strain interactions, etc., predicted cor- 
rectly? Or arc there compensating assumptions where one pro- 
cess is overpredicted at the expense of another and yet the end 
predicted result for the mean flow agrees with experiment? 
The lack of complete sets of data for higher moments, dissi- 
pation, and pres sure- velocity correlations for various flows 
has prevented detailed verification of closure models for the 
various processes that have to be modeled. 

The objective of this paper is to use the recently obtained 
and comprehensive experimental data of Shabbir and George 
(1987) and Shabbir (1987) on the axisymmetric buoyant plume 
to assess the various closure relations proposed for the kinetic- 
energy/ dissipation and the algebraic stress models for buoy- 
ancy-dominated flows. The usual approach is to solve the 
modeled differential equations numerically, and then compare 
the computations with the experiment. However, this method 
docs not help pinpoint the drawbacks in the various terms of 
the models. In this paper, instead of the usual approach, cor- 
relations obtained from measured velocity and temperature 
are used directly to verify the closure hypotheses for the tur- 
bulent transport of momentum, thermal energy, and turbulent 
kinetic energy. 

2 Experimental Data 

The data used were taken in an axisymmetric buoyant plume 
by Shabbir and George (1987) and Shabbir (1987), who meas- 


Contriliutcd by the Heat Transfer Division and presented at the National Heat 
Transfer Conference. Pittsburgh, Pennsylvania. August 9-12. 1987. Manuscript 
received by the Heat Transfer Division March 3. 1988; revision received No- 
vember 10, 1989. Keywords: Modeling and Scaling, Plumes, Turbulence. 
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ured velocity and temperature fields at several vertical levels 
above a heated source of air. Here wc briefly summarize their 
experimental technique and results. 

The three-wire probe used consisted of a cross-wire and a 
temperature wire. Thus, the instantaneous values of the two 
velocity components (vertical and radial) and temperature were 
measured. The axisymmetry of the flow was established by 
using an array of 16 thermocouples and also by rotating the 
cross-wire by 90 deg. Profiles for the correlations between the 
velocity components and velocity components with tempera- 
ture through the fourth order were determined from the in- 
stantaneous measurements. 

Source conditions were continuously monitored in order to 
calculate the rate at which buoyancy was added at the source. 
The source Grashof number was 5.5. By integrating the mean 
energy equation, an integral constraint can be obtained for a 
buoyant plume. For a neutral environment this constraint im- 
plies that the rate at which buoyancy crosses each horizontal 
section is constant and must equal the rate at which buoyancy 
is added at the source, i.e., the ratio 

JT = JF [ 2t j 0 S^(UAT+ut)rdi^ (1) 

must be unity (F 0 is the source buoyancy). This integral con- 
straint was satisfied within 7 percent. 

The correlation profiles at various heights were found to be 
similar in the coordinate rj = r/z {z accounted for the virtual 
origin) when the velocity is scaled by U s F\ n z~ ul and the 
temperature is scaled by T t = F^z" sn /g0. The measurements 
agreed well with the earlier study by George et al. (1977), who 
measured only the temperature and the vertical component of 
velocity. The scatter in the measurements of higher moments 
is typical for such flows and is also present in previous ex- 
periments, such as those of George et al. (1977). The primary 
reason for the scatter is that slow time scales of the flow require 
much longer averaging lime for the higher moments in order 
to obtain the same statistical convergence as for the mean 
quantities. Other errors in the measurements arise from the 
flow reversal on the hot wire — a phenomenon most likely to 
occur toward the outer edges of the flow where local turbulent 
intensities arc considerably higher. These arc discussed in Shab- 
bir and George (1987). 

The various correlations in similarity variables were fitted 
with curves using a least-squares fitting procedure. This rep- 
resentation allows easy evaluation of the terms in the governing 
equations and closure formulations when they are cast in sim- 
ilarity variables. Using these profiles the balances for the mean 
momentum and energy differential equations were carried out 
to check whether the flow satisfied the equations of motion it 
is supposed to represent. Within the thin shear layer and the 
Bousstnesq assumption the mean momentum and energy equa- 
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Fig. 1 (j) Balances of mean energy and momentum equations (taken 
from Shabbtr and George, 1987) 


lions can be respectively written as 

U-^~+ K ir == ~“7' (™>) -Z&& T ( 2 ) 

dz dr r dr 


3z dr 


--j{rvl) 
r dr 


( 3 ) 


Since all the quantities appearing in these equations are meas- 
ured, their profiles were substituted to see whether the meas- 
urements balance the equations. Figure l(o), taken from 
Shabbir and George (1987), shows that the experiment satisfies 
this nontrivial test within 10 percent. An error of such mag- 
nitude Is typical of turbulent shear flows. 

The dissipation of mechanical energy was determined by 
balancing the turbulent energy equation 



Fig. 1{b) Mechanical and thermal dissipation profiles! Full lines are 
experimental; the chained line Is from model equation (10), which was 
solved (or « with all other quantities taken from experiment 






+P+G-C 


( 4 ) 


where P = — TtjUj dU/dxj is the mechanical production and 
G = - Pgiu? is the production due to buoyancy. Each term 
except for the dissipation « and the pressure transport pUi is 
determined from the experimentally determined correlat ions . 
The pressure transport was evaluated from JnTj/p = 

5, a formula given by Lumley (1978). Although this closure 
relation has not been verified experimentally, it was felt that 
since the pressure transport is significant, some correction 
should be included rather than simply neglecting it, as is often 
done. The dissipation determined from the balance of the 
turbulent kinetic energy equation is shown in Fig. 1(a) as a 
solid line. 

By a similar procedure the dissipation of the mean-square 
temperature r 2 is determined from 






( 5 ) 


All terms are evaluated from experimental data and the re- 
sulting thermal dissipation is shown in Fig. 1(b). 

The time scales q 1 /* and f 2 /*, for the relaxation of the me- 
chanical and thermal dissipation, respectively, are shown in 
Fig. 2, along with their ratio 

R= (?/<,)/(?/<) ( 6 ) 


Nomenclature 


F m = buoyancy flux, equation (I) 
g = acceleration due to gravity 
G — turbulence production from 
buoyancy 

k = turbulent kinetic energy 
p = fluctuating pressure 
P = turbulence production by mean 
flow 

Pr r = turbulent Prandtl number 


r = radial coordinate 
R « time scale ratio, equation (4) 
t = fluctuating temperature 
T *= mean temperature 
u * fluctuating axial velocity com- 
ponent 

U — mean axial velocity component 
u = fluctuating radial velocity com- 
ponent 


z ~ vertical coordinate 
P = coefficient of thermal expan- 
sion 

< = dissipation of mechanical en 
ergy 

= dissipation of mean-square 
temperature 

- turbulent eddy viscosity 
P = mean density 
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Transactions of the ASME 







These time scales appear extensively throughout the model 
formulations and will be further discussed in the next sections. 


3 Assessment of Closure Hypotheses of k-< Model 

The form of the k-t model, which is considered to be the 
standard one, is that used by Launder and Spalding (1974). 
In this model the Reynolds stress is given by 



and the heat flux by 


p r ar 


where » T - and C , = 0.09 (see Launder and Spalding, 

1974). 


W - - <* r /Pr r ) — 

Taking Pr r = 1.0, the right-hand sides of the above equations 
were evaluated experimentally. These are compared with meas- 
ured values of uv, ul , and VI in Figs. 3-5. The points are 
experimental values and the chain lines are from the model. 

The modeled values of Vu and V7 compare reasonably with 
the experimental profiles except in the outer portion of the 
curves. On the other hand the modeled profile of vertical heat 
flux ul is much smaller than the experimental one. It is well 
known that the simple gradient models given by equations (7) 
and (8) with an isotropic eddy viscosity arc inadequate for 
determining streamwise turbulent momentum and heat fluxes. 
Usually these quantities do not influence the prediction for 
shear layers since only the radial fluxes are important in these 
flows. However, in the case of the buoyant plume the flux 
177 is a dominant production term in the turbulent kinetic energy 
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equation and its correct calculation is very important for ac- 
curate prediction of k. 

The diffusional transport in the kinetic energy equation (4) 
for a thin shear layer is modeled as 


1 — y l dk 

-pu=-r r —- 

2 p or 


(9) 


Using pv/p = -vqf/S from Lumley (1978) gives vq 1 / 2 = 
— {S/Z)v T Bk/Br from which the result, with the right side eval- 
uated from experimental results; is shown in Fig. 6. It is seen 
that the predicted and experimental data peak at different 
radial locations; however, the predicted magnitude is more 
accurate, which indicates that the pressure diffusion needed 
to be taken into account. 

In the k-t model the dissipation is calculated from 




( 10 ) 


where t'uj = -{vf/o^Bt/Bxj o t = 1.3, C ti = 1.44, and 
C< 3 = 1.92 as given by Launder and Spalding (1974). In order 
to get an indication of the validity of equation (10), it was 
numerically solved for the dissipation t with all other quantities 
needed to evaluate the coefficients determined from the ex- 
perimental correlations. The result is shown in Fig. 1(6} and 
it is seen that it compares reasonably well with the curve ob- 
tained from balancing the turbulent kinetic energy equation 
with experimental data. 

Launder et al. (1972) showed that the standard k-t model 
yields a solution for the axisymmetric jet that overpredicts the 
spreading rate by about 30 percent. The standard k-t model 
also does not correctly predict the axisymmetric buoyant plume 
(Hossain and Rodi, 1982). Proposals have been made (Pope, 
1978; Hanjalic and Launder, 1980) for modifying the dissi- 
pation equation, based on arguments concerned with vortex 
or eddy structures characteristic of axisymmetric flows. The 
modified equation produces more dissipation, thus decreasing 
the turbulent eddy viscosity, which results in a smaller spread- 
ing rate of the flow. Here we use the empirical correction given 
by Rodi (1972) where C« 2 = 1.92 (1-0.035 H) with H = 1 {y E / 
U m )d(J m /dx I 0 " 2 and where U m is the maximum velocity and 
y E is the distance from the centerline to the edge of the shear 
layer. This correction decreases the destruction term in the 
dissipation equation, hence producing an increased dissipation. 
However, when this correction is used in equation (10) there 
is very little change in the solution for « when experimental 
data are used for the other quantities in the equation. This is 
probably due to the approach taken here, which does not allow 
for the nonlinear interactions between the various terms in the 
closure. If the kinetic energy and dissipation equations are 
solved simultaneously, then the axisymmetric correction will 
produce a significant change in the solution of the k-t model. 


4 Assessment of Closure Hypotheses for Algebraic 
Stress Model 

Chen and Rodi (1975), Tamanini (1978), Chen and Chen 
(1979), and Hossain and Rodi (1982) have made predictions 
for the buoyant jet using algebraic stress models. Many of the 
ideas used in these models for calculating buoyant flows orig- 
inated with Launder (1975. 1978). Algebraic stress models are 
obtained by simplifying the convective transport equations for 
Reynolds stresses and heat fluxes so they are no longer dif- 
ferential equations. The dynamic equation for the Reynolds 
stress tensor is 



where P~ - -u,-u k BU/Bx k - u/i k BU/Bx k is the mechanical 
production and G<, = uj - figj is the buoyancy 

production. The left side represents convection minus diffu- 
sional transport, the dissipation is assumed isotropic, and the 
last three lines rep resent the closure formulation for 
p(du/dxj+ duj/dx d )/p given by Launder ct al. (1975) and 
Launder (1975, 1978). Launder assumes: (I) an equilibrium 
situation where convection is balanced by diffusion (C - D 

- 0) and production is balanced by dissipation (/> + G — t 

— 0); (2) the second and third terms (third line) in the rapid 
part of the pressure-velocity correlation are negligible; the 
coefficient C 2 is adjusted so that the first term approximates 
the entire rapid pan; (3) the parameter C 3 is taken equal to 
C 2 . After applying ali the assumptions 


2 C, + C 2 - 1 

3 C, 


*5 V - 


1-C 2 * 

c, t 


(P^G^) 


( 12 ) 


where c, = 2.2 and C 2 = 0.6. It should be pointed out that 
in free shear flows the equilibrium condition (C — D ** 0 and 
P + G - t - 0) only applies in the outer portion of the flow. 
Also, Zeman and Lumley (1976) found C 3 = 0.3, after applying 
ali the constraints applicable to determining the contribution 
of buoyancy to the pressure-strain correlation. 

The dynamic equation for the heat flux is 




+ +C 3< 0*/ (13) 

where the first line on the ri ght side is the production and the 
second line is the closure for pBt/BxJp. Neglecting convection 
and diffusion (C — D = 0) and the third term in the second 
line, equation (13) becomes 

■ 1 k f BT BUi 

5?= r u * ^ - u- c »w*?\ 

(14) 

where C„ = 3.0, C* = 0.5, and Cj, = 0.5. Zeman and Lumley 
(1976) show that C = 0.8 and C 3/ = 0.2 from theoretical 
considerations. 

Neglecting the convection and diffusional transport in equa- 
tion (5) and eliminating *, with equation (6) gives 


< Bx, 


OS) 
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which was given by Launder (1975, 1978) and used by Hossain 
and Rodi (1982). 

Chen and Rodi (1975) and Chen and Chen <1979) used the 
differential cquatjpn (5), with eliminated by using equation 
(6) to determine t 2 rather than using equation (15). In cither 
case R is taken to be a constant equal to 0.8 (Hossain and 
Rodi, 1982; Chen and Rodi, 1975; Chen and Chen, 1979; 
Launder, 1975, 1978). It is seen in Fig. 2 that the experimentally 
determined value of R is much lower with an average value 
across the profile of roughly 0.25. Launder (1978) cites ex- 
perimental evidence for R being in the range of 0.5 to 0.8. 
However, he found that the algebraic stress relations agreed 
best with an experiment for a stably stratified homogeneous 
shear flow with R = 0.8. Hence, that value has been adopted 
in the algebraic stress models. The experimental results of 
Shabbir and George (1987) indicate/? is much lower for strongly 
buoyant flows. When the algebraic stress model is applied to 
this experiment with R — 0.8, the results arc very poor for 
r 2 . Therefore, in the following evaluation of the algebraic stress 
model, the experimentally determined profile for R (Fig. 2) is 


used. 

Equations (10), (12), and (13) represent a system of algebraic 
equations that caii be solved for uiu j% uj % and r 2 . Employing 
the thin shear layer approximation, where only gradients in 
the radial direction are retained, Hossain and Rodi (1982) give 


2C, + C 2 -1 

3 Cj 


k + 




(16) 



Fig. 7 Mean square temperature 



r/z 

Fig. 8 Radial Reynolds stress tr* 


. i-C 2 */ -jdU _Y 


75_ 2 Q + Cj- 1 
3 C, 


(17) 

(18) 


57= '2 — f — uv — ( 1 — Ci,)57 — + ( 1 — CxIPgt 1 

C u ( L dr- . dr J 


_ -.1 k-BT 

?= - 277 - 57 ^ 

. < -dr 

from which v, = C r k?/t where 
- _ 2 (t-q,)(C, + C 2 -l) 

3 <3 


( 


, 1 k - 3 773 A 

+ C u t gP dU/dr) 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 


The system, of equations (16)-{21) was solved to determine the 
Reynolds stress and heat flux components with (7, T, k t e, and 
R given by the experiment. The following values for the con- 
stants were used: 

C,=2^, C 2 = 0.6, C 3 * 0.6 

C w -3.0, C*=0.5, Q,= 0.5 

The value of C^ t which appears in the eddy viscosity relation 
r, = CJP/t and is given by equation (22), is roughly 0.125 
and is reasonably constant across the flow. This value is con- 
siderably larger than the value of C* « 0.09 in the standard 
k~t model. Rodi (1972), to correct for the discrepancies in the 
prediction for the axisymmetric jet, developed an empirical 
correction to The parameter C # is replaced by (1-0 .465//) C,, 
where H = I {y^/U^dU^dxi U m is the maximum velocity, 
and y £ is the distance from the centerline of the edge of the 
jet. Hossain and Rodi (1982), Chen and Rodi (1975), and Chen 
and Chen (1979) used this correction in their predictions for 
turbulent buoyant jets. When the correction is applied, we get 


approximately 0.09, which is the value of for the standard 
k-( model. 

The question we are asking is, “given the turbulent energy, 
dissipation, velocity, and temperature, does the proposed al- 
gebraic stress expression correctly predict the Reynolds stress 
and heat flux components?** Figure 8 shows the radial Rey- 
nolds stress determined from equation (18). It is seen that the 
predicted value x^/k = 0.53 is a little smaller than the exper- 
imental curve in the center portion of the plume, but agrees 
quite well with experiment in the outer portion. The predicted 
shear stress uu, the axial heat flux i/7, the radial heat flux u7, 
and the mean squared temperature r 2 are shown in Figs. 3, 4, 
5, and. 7, respectively. Again the points are experimental data 
and the broken lines arc from the model. It is seen that the 
shear stress uv and radial heat flux vl are predicted reasonably. 
However, the vertical heat flux ul and temperature fluctuations 
* r 2 arc predicted poorly and have incorrect shapes; unlike the 
experimental values they go to zero near the origin. 

Equation (21) gives r 2 proportional to the radial temperature 
, gradient, which is zero at the centerline. Then since r *= 0 at 
r = t), equation (19) gives a7 - 0 at r - 0. In order to obtain 
nonzero values; for ul and t 2 at. the centerline from the model 
equations (12), (14), and (15), terms containing the axial gra- 
dient, i.e., dU/dz and dT/dz, were retained. These terms were 
added to equations (19) and (21) and the system of equations 
was- solved again. Although the centerline "values of ut and 
r 2 were found to be nonzero, the predictions still d<kxea sed to 
relatively small values. near the centerline. 

Another possibility for this behavior is the neglect of ad- 
vection and diffusion terms in the model. Gibson and Launder 
(1976) have proposed the following model for these terms: 


{.C-D)zpj= ^ < P+G-t ) 


(23) 


57 57 

where P t is the production term in the t 2 equation. These were 


(C-0)^= (/>,-*,)+ ~ (F+G-e) (24) 
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Fig. 9 Transport of temperature variance t* 


incorporated in equations (I6)-(21) and the resulting set of 
nonlinear coupled algebraic equations was solved simultane- 
ously. The results did improve the prediction for the vertical 
heat flux ul and temperature variance r 2 but the comparison 
for radial heat flux and shear stress became worse. As noted 
by Gibson and Launder (1976), the above model is not good 
near an axis of symmetry. This is why, by incorporating them 
in the original model, no overall improvement in the prediction 
is achieved. 

Ghen and Rodi (1975), Tamanini (1978), and Chen and Chen 
(1979) use {he differential convective-transport equation (5) to 
determine t 1 in their predictions of buoyant jets. Equation (6) 

was used to eliminate Thus, the final form of the r 2 equation 
becomes 


vf+vf 



(25) 


This equation was numerically solved for the temperature var- 
iance r 2 with all other quantities needed to evaluate the cocf- 
'ficiehts deteimined from experiments; The value of C f was 
taken as 0.13. The best agreement, as shown in Fig. 7, was 
achieved with R = 0.35. When the average experimental value 
of R « 0.25 is used, the prediction peaks at about 6.0 (tj = 
0.04) as compared to the experimental Value of? of about 8.0 
(ij ~ 0.04). When the standard value of R = 0.8 is used, ri 
is overpredicted by a factor of four. 

When f 2 is calculated from the convective-transport equation 
(5), the diffusive transport is given by the simple gradient 
closure • - * 



' « dr 


(26) 


wth C, «= 0J3 as given by Chen and Rodi (1980). The pre- 
diction for trf 2 , using Experimental information to evaluate the 
right-hand side of equation (26), is shown in Fig. 9. The pre- 
dicted curve peaks somewhat above and toward the centerline 
as compared to the data. 

Finally, we ask that if the models do not depict the axial 
heat flux 177 and the temperature variance ? correctly, then 
why do the predictions such as made by Hossain and Rodi 
(1982), Chen and Rodi (1975), Tamanini (1978), and Chen and 


Chen (1978) show reasonable agreement with the experiment 
for the mean velocity and buoyancy? The answer to this is 
that with R - 0.8 the temperature variance r 2 from equation 
(21) or (25) is too large. This makes the vertical heat flux ul 
from equation (19) large enough so that the mean velocity and 
buoyancy are reasonably predicted. 

5 Summary and Conclusions 

The experimental data on buoyant plumes were used to 

evaluate various closure relations for turbulence transport. The 

objective was not to propose new models, but to evaluate the 
closure schemes proposed by other workers for buoyancy- 
dominated flows. The closures evaluated were those used in 
the k-t and algebraic stress models. The results are summarized 
below. 

1 The closure relations of the Jr-< model compare reason- 
ably with experimental data, except for the axial turbulent 
transport, which is drastically underpredictcd. The axial heat 
flux governs the production due to buoyancy in the kinetic 
energy and dissipation equations and its correct prediction is 
very important. This is a probable reason why the results of 
Hossain and Rodi (1982) from the k-t model under predict the 
spreading rate for the plume by 10 percent even when axisym- 
metric jet corrections are included. 

2 The ratio R of the time scales, which is used to determine 
the dissipation of the mean squared temperature in the alge- 
braic stress model, was found to be considerably different from 
the accepted value of R - 0.8. Apparently R is not a universal 
constant, but can vary from flow to flow and is influenced by 
the strength of the buoyancy present. From the experimental 
data on a plume it appears that R - 0.25 for strongly buoyant 
flows. 

3 The closure equations for the shear stress and radial heat 
flux of the algebraic stress models also compared well with 
experiment but arc not better than the simple gradient closures 
used in the k-t model. The axial heat flux and mean squared 
temperature are predicted poorly in the central core of the flow 
and had incorrect trends. This drawback could be attributed 
to the assumption of local equilibrium, which resulted in the 
neglect of convection and diffusion terms in the transport 
equations for Reynolds stress and heat flux. However, no 
substantial improvement was.achieved by keeping the second- 
ary derivatives or by incorporating the model for the convec- 
tion and diffusion terms. Therefore, the full dynamic equations 
for Reynolds stress and heat flux with convection and diffusion 
are required to predict the axial heat flux and temperature 
variance properly. 
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